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ABSTRACT 

We examine deep XMM-Newton Reflection Grating Spectrometer (RGS) observations of the 
X-ray luminous galaxy cluster A1835. For the first time in a galaxy cluster we place direct 
limits on turbulent broadening of the emission lines. This is possible because the coolest X-ray 
emitting gas in the cluster, which is responsible for the lines, occupies a small region within 
the core. The most conservative determination of the 90 per cent upper limit on line-of-sight, 
non-thermal, velocity broadening is 274 km s~ 1 , measured from the emission lines originating 
within 30 kpc radius. The ratio of turbulent to thermal energy density in the core is therefore 
less than 13 per cent. There are no emission lines in the spectrum showing evidence for gas 
below ^3.5 keV. We examine the quantity of gas as a function of temperature and place a 
limit of 140 M yr~ 1 (90 per cent) for gas cooling radiatively below 3.85 keV. 

Key words: X-rays: galaxies — intergalactic medium — galaxies: clusters: individual: Abell 
1835 — cooling flows 



1 INTRODUCTION 

Little direct information is known about the dynamical state or 
microphysics of the intracluster medium (ICM). Simulations of 
galaxy clusters within a cosmological environment predict that as 
matter accretes along filaments, or as clusters merge, turbulent mo- 
tions should be set up in the ICM. This energy should then cascade 
from large to small scales where it can be dissipated. Simulations 
have predicted that even in the most relaxed clusters there are sub- 
stantial flows of ICM (e.g. Evrard 1990). In these simulations the 
fraction of pressure support in gas motions is typically 5 to 15 per 
cent. Most simulations, however, have not included physical pro- 
cesses which may be relevant, such as viscosity or magnetic fields. 

In addition, it is observed that the active galactic nuclei (AGN) 
in cluster cores inject jets and bubbles of relativistic plasma into 
their surroundings. The feedback of energy into the cluster is 
thought to combat the high radiative losses of energy which would 
otherwise lead to substantial cooling rates (for reviews see |Peter^| 
son & Fabian 2006 and McNamara & Nulsen 2007 ). It is predicted 
that feedback could induce gas motions i n the core in the range 
500- 1000 kms -1 i Bruggen et al.||2005 i. Some models also re- 
quire turbulence in order to heat the ICM using the energy in in- 
flated cavities (e.g. |Scannapieco & Brtiggen 2008). Strong turbu- 
lent motions would provide significant non-thermal pressure sup- 
port. This would bias the determination of cluster mass profiles 
measured under the assumption of hydrostatic equilibrium. The 
measurement of turbulence is thus important for mass determina- 



tion, as well as AGN feedback. 

Evidence for resonance scattering has also been used to place 
upper limits on turbulence in elliptical galaxies. In NGC 4636, |Xu] 
J2002| examined the strength of the resonantly scattered 15/ 



et al. 



Fe XVII line relative to the line at 17.1 A, inferring scattering. As 
turbulence broadens emission lines and prevents scattering, they 
obtained an upper limit on the turbulent velocity dispersion of 10 
per cent of the sound speed. For the same object, Wern er"et al.| 
( 2009 1 calculated a limit of 25 per cent of the sound speed, with an 
implied maximum of 5 per cent of energy in turbulent motions. 

The Perseus cluster is an interesting test case. The tens of 
kpc long quasi-linear filaments ( Fabian et al. 2003b I imply laminar 
flow, significant viscosity in the ICM and low levels of turbulence 
on the scale of the filaments. However, the lack of evidence for 
resonant scattering in the strongest 6.7 keV line of He-like iron in 
Perseus argues for velocity spread of at least half the sound velocity 
(Chur azov et al.|2 004). Reb usco et al.|p005] l also require motions 
of 300 km s~ 1 on scales of 20 kpc to account for the overall metal- 
licity profile in Perseus. 

Comparing the gravitational mass profiles derived from opti- 
cal and X-ray data in M87 and NGC 1399, |Churazov et al.|p07)8) 
obtained combined nonthermal pressures of ~ 10 per cent of the 
thermal gas pressure. Reynolds et al. (2005) simulated the evolu- 
tion of cavities in the ICM, requiring viscous processes in order 
for the cavities to retain the observed spherical cavity shapes and 
prevent them from being shredded by instabilities (Reynolds et al. 
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Table 1. Single temperature RGS spectral fitting results. The results with 
RGSXSRC use that model to account for the spatial extent of the source. 



Parameter 



Without RGSXSRC With RGSXSRC 



kT (keV) 

O(Z ) 

Ne(Z a ) 

Mg(Z ) 

Si(Zc) 

Fe(Z ) 

normalization (cm - 
C statistic 



3.67±0.16 
0.17±0.03 
0.25 ±0.04 
0.17±0.08 
0.19±0.06 
0.18±0.02 
0.0104 ±0.0001 
3947.9 



3.71 ±0.16 
0.21 ±0.03 
0.30 ±0.05 
0.20 ±0.09 
0.19 ±0.06 
0.20 ±0.02 
0.0107 ± 0.0001 
3956.4 



|2005l|Dong & Stone|2009) . A lower limit of 10 per cent of the total 
pressure in turbulence was found in the Coma cluster by Sc huecker| 
|et al.| ([2004). The Coma cluster is unrelaxed with two central galax- 
ies. It is unlikely that it can be directly compared to relaxed clusters 
or ellipticals. In summary, the evidence for turbulence in relaxed 
clusters is uncertain. 

A 1835 (at z = 0.2523) is the most luminous cluster in the BCS 
sample (Ebeling et al. 1998), with an inf erred cooling rate of 



1000 Mr.-, yr in the absence of heating (Allen et al. 



1996). The 



optical-spectroscopic-der ived rate of star form ation in the central 



galaxy is 40-70 M s yr 1 1 Crawford et al. 1999 1. The IR luminosity 



of the object is ~ 7 x 10 Lq i Egami et al.|2006 i, implying a star 



formation rate of ~ 125 Mq yr . In addition there is ~ 10 M : 
of molecular gas in the cluster (Edge 2001 ). 

Despite the evidence for star formation, cool gas and the high 
X-ray luminosity, little X-ray emitting gas was seen in this cluster 
below 1-2 k eV when it was observ ed using the XMM-Newton RGS 
instruments (Peterson et al. 1|200l| . Less than 200 M@ yr 1 cools 
below 2.7 keV (90 per cent confidence). AGN feedback, at a level 
of 1.6 x 10 45 erg s -1 , prevents most of the X-ray gas from cool- 
ing. Later observations confirmed the lack of X-ray cool gas is a 
common feature of clusters (e.g. Peterson et al. 2003}. 

In this paper we examine deep RGS observations of A 1835 
totalling 254 ks, placing limits on velocities in the gas and the 
amount of cool gas detectable in X-rays. We use a Hubble c onstant 
of 70kms~' Mpc~' and the relative Solar metallicities of Anders 
|& Grevesse| ( [T989l >. 1 arcsec corresponds to 3.94 kpc at z = 0.2523 
Errors quoted are 1<T unless stated otherwise. 



2 RGS SPECTRA 

We examined three XMM-Newton RGS observations (0098010101, 
0551830101 and 0551830201). We extracted the RGS spectrum 
from the datasets using the standard SAS 9.0 RGSPROC tool. We 
used a 90 per cent PSF width extraction region (corresponding 
to a half-width of approximately 30 arcsec), and a 90 per cent 
pulse-height distribution energy selection. To eliminate flares, we 
removed time periods with a count rate greater than 1 s on CCD 
9, selecting events with FLAG set to 8 or 16 and a cross dispersion 
angle greater than 1.5 x 10~ 4 . 

Background spectra were generated from the observations, ex- 
tracted from the region beyond a 98 per cent of PSF width. We 
tested that the backgrounds were essentially source free by com- 
parison with template background spectra. The wavelength binning 
option was used so that we could combine the spectra from both 
RGS instruments using RGSCOMBINE. We combined the spectra, 
backgrounds and responses from all observations, treating the first 
and second order spectra separately. 



2.1 Single temperature model 

We show in Fig. [T] the first order combined spectra with the best 
fitting thermal model. We used the SPEX version 2.00. 1 1 thermal 
model (Kaastra & Mewe 2000), converted to an APEC format ta- 
ble for XSPEC, with Galactic absorption fixed at 2.04 x 10 20 cm~ 2 
(Kalb erla et al.||2~005| l. We fitted both spectral orders simultane- 
ously, using a wavelength range of 7 to 28A and 7 to 16A for the 
1st and 2nd orders, respectively. The O, Ne, Mg, Si and Fe metal- 
licities were allowed to be free in the fit with the other elemental 
abundances tied to Fe, except He which was fixed to Solar. Ni and 
Ca lines in the fitted spectral region are too weak to fit separately. 
We allow a normalization offset between the two orders (the 2nd 
order gives a best fitting normalization 3 per cent greater than the 
1st order spectrum). We fit the spectrum minimizing the C statis- 
tic on the ungrouped spectra in XSPEC 12.5.1 ( Arnaud 1996}. The 
thermal broadening option was enabled when fitting the spectra. 

The C statistic does not provide a goodness of fit, but the 
XSPEC goodness command calculated 50.9 per cent of Monte Carlo 
realisations of the best fitting model have a better fit statistic than 
the real data, indicating a good fit. If the spectra are fit with the same 
model but using % 2 statistics, the reduced % 2 is 1-02 = 2253/2218 
(grouping to have at least 20 counts per channel). Shown in the 
lower panel of Fig. [T] are the rebinned residuals of the single tem- 
perature model to the first and second order spectra. The best fitting 
model parameters are shown in Table[T] (under the column without 
RGSXSRC). The low Fe abundance is interesting. This could be sig- 
nificantly affected by how much continuum emission is from hotter 
plasma (> 5 keV) outside the cluster core, which is not well con- 
strained by lines in the RGS band. We note that the line complexes 
atll.4andl2.2A are underpredicted by our model, which we could 
not match by adding additional temperature components. 

Fig. [T] shows the positions of emission lines that are seen by 
eye in the spectrum and the positions of lines that are not seen. All 
ionization states of Fe seen are Fe XXIII and above. The detections 
of Fe XXIII are not definite as they are blended with other lines 
(Fe XXIV and Ne X). If we plot the total line emissivity of ioniza- 
tion states seen and not observed in the spectrum as a function of 
temperature (Fig. [2}, there are no obvious lines that indicate gas at 
temperatures below ~ 3 keV. This is also seen by comparison with 
the isothermal spectral models in the middle panel of Fig. [T] 

If we add a further component with the same metallicities as 
the first, it has a best fitting temperature of 0.58^jJ'jg keV and an 
XSPEC normalization of 6.1^2g x 10~ 5 cm~ 5 , but is only signifi- 
cant at the 90 per cent confidence interval. 



2.2 Smoothing model (RGSXSRC) 

As the cluster is not a point source and the RGS instruments are slit- 
less spectrometers, the observed spectra are broadened according to 
the spatial extent of th e source. This broade ning is approximately 
AX = (0.124/m)Ae A |Brinkman et al.|l998| , where m is the spec- 
tral order and A6 is the half energy width of the source in arcmin. 
The effect is relatively small in A1835, relative to nearby clusters. 

Broadening was modelled using the RGSXSRC model in 
XSPEC with an input Chandra image in the 0.8 to 1 keV band. 
The position angle of XMM was similar between the observations 
(~ 288°). RGSXSRC assumes that the source has the same spectrum 
as a function of position, so is not completely correct if different 
emission lines come from different radii in the cluster. 

The second column in Table [T] shows the small effect of this 
smoothing on the fit parameters. The fit quality is slightly poorer 
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Figure 1. (Top panel). Combined RGS 1 and 2 spectrum from the 90 per cent extraction region and best fitting single-temperature model. The data have been 
rebinned to have a signal to noise ratio of 10 and divided by the effective area as function of wavelength. The spectral feature close to 16A seen in the spectral 
model next to Fe XVIII is O VIII. (Centre panel) Isothermal Solar metallicity spectral models at different temperatures smoothed with the RGS response. 
(Bottom panel) Rebinned residuals of the model fit to the data. The magenta points at short wavelength are the residuals from the second order spectrum. 



after accounting for smoothing, suggesting that the emission lines 
may be emitted on a smaller scale than the broader X-ray emission 
(Section]?}. The goodness command gave 54.6 per cent of realisa- 
tions of the best fitting model with a better fit than the data. 



2.3 Limits on velocity broadening of the spectra 

Since the lines are remarkably narrow, and the source extent is 
small, we can limit the velocity broadening. We fitted the spectra 
in Section[2~T]with the BVAPEC model in XSPEC assuming that the 
spectra were broadened by the thermal motion of the ions in the gas 
(version 12.5.0ah of XSPEC fixed an error in the line widths) and 
line-of-sight velocity broadening added in quadrature. The SPEX 
spectral lines were used in the BVAPEC model. 

We show the change in fit statistic as a function of velocity in 
Fig. [3] If we conservatively treat the cluster as a point source and do 
not use the RGSXSRC model we obtain a 90 per cent upper limit of 
274 km s _ 1 . This limit is improved to 1 82 km s _ 1 with the addition 
of the broadening model. Examining just the Fe-L spectral region 
between 9.6 and 12.8A (rest), we obtain a limit of 214kms~' , and 
for the 1st order spectrum between 13.6 and 22.4A, containing the 
strong O VIII line, a limit of 380 km s~' . We confirmed that such 
constraints are readily achievable with simulated spectra. 



2.4 Limits on cooling gas 



We can limit how much gas can be cooling within temperature 
bins. We constructed a model using six cooling flow model VM- 
CFLOW components (constructed using SPEX models) with temper- 
ature ranges of 10 -> 5.6^ 2.8 -»-1.4-»0.7-> 0.35 -> 0.0808 keV. 
We fixed the components to have the same metallicities, allowing 
individual elements to vary as in Section [2~T| The maximum cool- 
ing rates in M Q yr~ 1 , assuming isobaric cooling, were free param- 
eters in the fit. The results are shown in Fig. [4] with and without the 
RGSXSRC smoothing component. 



For comparison, we also fitted a model made up of an isother- 
mal component plus a model cooling from its temperature to zero. 
The best fitting temperature of the isothermal component was 
3.85 keV (close to the single temperature fits). We obtain cooling 
rates from this temperature to zero of 70+^ M e yr" 1 . This rate 
is shown by the shaded bar in Fig. [4] The 90 per cent upper limit 
is 140 Mq yr~' . A Markov Chain Monte Carlo analysis produces 
very similar values. Comparing an isothermal to isothermal plus 
cooling flow model with an F-test and x 2 statistics gives a chance 
improvement of the x 2 of 17 per cent. 



© 0000 RAS, MNRAS 000, 000-000 



L4 J.S. Sanders et al 



10-' 4 F 




Temperature (keV) 

Figure 2. Total emissivity of spectral lines from different species as a func- 
tion of temperature. The solid lines are observed in the RGS spectrum of 
A 1835 and the dotted lines are not. 

3 CLUSTER PROPERTY PROFILES 

For comparison with the RGS results we show projected Chandra 
and XMM EPIC-MOS temperature profiles and deprojected Chan- 
dra temperature, density and cooling time profiles. (Fig. [5}. We 
examined Chandra datasets 6680, 6881 and 7370 (total 170 ks, 
cleaned), reprocessed applying VFAINT filtering. Blank-sky back- 
ground spectra were used, scaled to the observed 9-12 keV count 
rates, providing a good match to the data in the field outskirts. EM- 
CHAIN was used to process each MOS dataset, removing bad time 
periods with a 10-15 keV rate > 0.35 s . Scaled blank-sky spec- 
tra were used as backgrounds. There was a small mismatch between 
the background and data in the outskirts which we modelled as a 
uniform model component. 

There is a good match between the XMM and Chandra pro- 
jected profiles, except where the PSF of XMM is important in the 
centre. In the core of the cluster the temperature drops to less than 
3 keV and the cooling time is around 5 x 10 8 yr. The entropy of 
the gas in the central region drops to 13keVcm 2 . The metallicity 
rises smoothly inwards from 0.2 to 0.4 Z . The temperature from 
the RGS spectral fit, 3.7 keV, is emitted from within 30 kpc radius 
(8 arcsec), essentially a point source for the RGS. 



4 DISCUSSION 

The limits we place on the broadening of the emission lines by 
non-thermal motions of the ICM are the first direct limits on turbu- 
lence in cool core clusters. As the cooler gas below 4 keV is located 
within 8 arcsec of the core (Fig. [5}, the conservative point-source 
limit of a line-of-sight velocity of 274kms~ 1 is more appropriate 
than the smaller RGSXSRC value. This velocity is 28 per cent of the 
sound speed at the measured temperature. 

Using equation 1 1 in Werner et al. (2009|l we can calculate the 
ratio of turbulent to thermal energy density, VjL [lm„jkl ', where 
Vi os is our measured line-of-sight velocity, fl is the mean particle 
mass, m p is the proton mass and kT is the temperature. We obtain 
a 90 per cent confidence upper limit of 13 per cent on the energy 
density in turbulence. 

Our upper limit is compatible with the values of 5 to 15 per 
cent of pressure in gas motions, calculated using high resolution 
Eulerian simulations ( |Lau et al.||2009| |Vazza et al~||2009| >. If vis- 
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Figure 3. Limits on velocity broadening of RGS spectra. Shown are the 
increases in fit statistic as the velocity is stepped from to 1000 km s , for 
the data fitted with the RGSXSRC smoothing model, without the model, and 
then fitting the 1 st and 2nd order spectra separately. 

cosity is significant in clusters, as indicated by the quasi-linear fila- 
ments in Perseus (Fabian et al. 2003b I, the real value is likely to be 
lower than simulations predict. The low turbulent pressure limit we 
measure from the core of this cluster is encouraging for the mea- 
surement of cosmological parameters using cluster mass profiles 
(e.g. |Allen et al.||2008| . Low turbulence implies that viscosity is 
significant. Sound waves generated by the AGN are then damped, 
heating the inner ICM ( [Fabian et al.|2003a| >. 

Our results also suggest that the feedback taking place in this 
cluster is not inducing large amounts of turbulence in the gas in 
the cluster core. The feedback must be gentle. The limit on turbu- 
lent velocities we obtain is much lower than the turbulence required 
in Perseus to prevent resonance scattering (Ch urazov et al.|2004[ >. 
However, the region studied in Perseus was large and may contain 
differential laminar motions. Shearing motions, rather than turbu- 
lence, may limit resonant scattering. The velocity limit we find is 
consistent with the velocity shear of 250 km s seen in Ha emis- 
sion for the central 3 1 kpc diameter nebula around the central clus- 
ter galaxy (Wilman et al. 2 006[ >. 

The RGS instruments have been used to measure line widths 
in extended Galactic objects (e.g. Rasmussen et al. 2001). We have 
shown here that the instruments are also capable of usefully con- 
straining turbulent velocities for X-ray bright and peaked clusters, 
of which this is the first. Non-cool-core clusters and mergers are 
not so spatially compact, making limits more difficult to measure. 
ASTRO-H will measure the ICM velocity distribution in a wider 
range of objects and with much greater precision. 1XO will map the 
velocity distribution over cluster cores and to higher redshifts. 

We are able to confirm the extremely low rates of gas cooling 
to low temperatures first measured by Pet erson et al.| ( [2001^ . The 
distribution does not, however, appear to be a powerlaw as seen in 
other objects (Peterson et al. 20031, at least below 4 keV. There 
is weak evidence for an additional residual cooling component 
at around 100M c;) yr _1 at 0.6 keV. This rate of cooling, and our 
other limits, are compatible with an o ptical-spectroscopic-d erived 
star formation rate of 40-70 Mq yr' '^Crawford et al.|| 1999> and 
of ~ 100 M yr _1 from UV fluxes ( jHicks & Mushotzky|2005), if 
the star formation is fed by cooling ICM. McNa mara et al.| ( 2006) 
suggest that the central galaxy hosts an extremely radiatively inef- 
ficient black hole the jet power of which is large enough to offset a 
large fraction of the cooling in this object. 
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Figure 4. Mass cooling rates as a function of temperature derived from 
RGS spectra. The solid points show maximum cooling rates within fixed 
temperature bins. The shaded region shows the rate and uncertainty using a 
model made up of an isothermal plus cooling flow component. 

Feedback is not completely balanced in this cluster, as is evi- 
dent by the high star formation rate, high IR luminosity and large 
amounts of molecular gas. A residual level of cooling of 10 per cent 
of the X-ray luminosity-inferred rate remains. Future detailed sim- 
ulations and modelling may allow us to place detailed constraints 
on the temperature, velocity and spatial distribution using distant 
luminous clusters like A1835. 



5 CONCLUSIONS 

For the first time we measure meaningful constraints of the ICM 
turbulent velocity from the width of the X-ray emission lines. We 
place a limit of 13 per cent on the turbulent energy density as a 
fraction of the thermal energy density. We also constrain the cool- 
ing rate of X-ray emitting gas in A1835 to a value comparable to 
the star formation rate. 
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